Magnetization- Switched Metal-Insulator Transition in a (Ga,Mn)As Tunnel Device 
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We observe the occurrence of an Efros-Shklovskii gap in (Ga,Mn)As based tunnel junctions. The 
occurrence of the gap is controlled by the extent of the hole wave-function on the Mn acceptor 
atoms. Using k • p-type calculations we show that this extent depends crucially on the direction of 
the magnetization in the (Ga,Mn)As (which has two almost equivalent easy axes). This implies one 
can reversibly tune the system into the insulating or metallic state by changing the magnetization. 

PACS numbers: 71.30.+h, 75.30.Hx, 75.50.Pp 



A very direct way to observe the Efros-Shklovskii 
(ES) gap, the soft gap induced by Coulomb correla- 
tions near the Fermi level of a Mott insulator 0, 0], 
is by means of tunnel spectroscopy. Such experiments 
were, e.g., performed on the (three-dimensional) non- 
metallic doped semiconductor Si:B [jj and on thin 
(two-dimensional) Be films0. While both of these ex- 
periments employed large area tunnel junctions and a 
metallic counter electrode, a more recent study employed 
Ge:As break junctions^. This latter approach avoids 
possible screening of the Coulomb correlations, but the 
mesoscopic character of the contact may complicate ex- 
traction of bulk Coulomb gap behaviour^- 

We have recently investigated the physics of a novel 
type of magnetoresistance, dubbed TAMR (tunneling 
anisotropic magnetoresistance) 0, E|. TAMR results 
from the dependence of the density of states (DOS) 
in strongly spin-orbit coupled ferromagnetic semicon- 
ductors, such as (Ga,Mn)As, on the direction of the 
magnetization of the material. In we reported a 
drastic (> 10 4 ) increase of the spin- valve signal in a 
(Ga,Mn)As/GaAs/(Ga,Mn)As tunnel structure on low- 
ering the sample temperature from 4.2 to 1.7 K, and spec- 
ulated that this behaviour might result from the opening 
of an ES gap. Here, we provide evidence that the high 
resistance state of the sample indeed corresponds to a 
soft-gapped Mott insulator. In these samples, the metal- 
to-insulator transition (MIT) is driven by a large varia- 
tion of the Bohr radius of a hole bound to a Mn-impurity 
when the magnetization of the layer is switched from one 
easy axis to the other. This assignment is supported by 
a k • p-type calculation of a hydrogen-like impurity in a 
ferromagnetic GaAs host, extending the successful mean 
field model for (Ga,Mn)As [H03- 

Our (Ga,Mn)As tunnel structure is shown in Fig. \%)p. 
From bottom to top, the Gao.94Mno.06 As (100 nm)/ 
GaAs (2 nm) /Gao.g4Mno.o6As (10 nm) trilayer stack has 
been grown by low temperature molecular beam epitaxy 
(LT-MBE) on a semi-insulating GaAs substrate and a 
120 nm undoped GaAs buffer layer. Both (GaMn)As 
layers are ferromagnetic with an as-grown Curie temper- 



ature of ~ 65 K and highly p-type due to the intrinsic 
doping arising from the Mn atoms. 

As seen in the optical micrograph of Fig. Q^, the layer 
stack is patterned into a square mesa of 100x100 fim 2 
by positive optical lithography, metal evaporation, lift- 
off and wet etching. The top contact is in-situ Ti/Au. 
Contact to the lower (Ga,Mn)As layer is established by 
a W/Au ring around the mesa, fabricated by metal de- 
position and optical lithography . 
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FIG. 1: a) Optical micrograph of the sample, b) Sample layer 
structure, c) schematic band diagram illustrating the barrier 
for hole transport, the depletion region near the interfaces, 
and the effective injector layers (i). d) 4.2 K and 1.7 K mag- 
netoresistance in zero field at 1 mV bias after preparing the 
sample by magnetizing in various field directions. 

The structure is intended for vertical tunneling trans- 
port experiments. A key element in this device design is 
that the LT-GaAs tunnel barrier contains many defects, 
and thus acts as an effective carrier sink, resulting in a 
thin depleted (Ga,Mn) As region near the barrier and thus 
band bending near the interfaces, as sketched in Fig. [Tfc. 
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We measure the resistance of the device depending on 
its magnetic history. Using a magnetocryostat equipped 
with three pairs of Helmholtz coils, we magnetize the 
sample along a specific angle (p in the sample plane. One 
of the magnetic easy axes of the (Ga,Mn)As layers, de- 
fined as <p = 0°, is along the [100] crystal direction. 
The second, perpendicular, easy axis is along [010] cor- 
responding to ip = 90°. As the field is swept to zero, the 
magnetization relaxes to the nearest easy axis, and this 
magnetization state determines the resistance. Fig. 
illustrates the results of such measurement [To| for many 
preparation angles <p. The 4.2 K-curve can be understood 
as resulting from TAMR 8]: the anisotropy in the mo- 
mentum dependent DOS with respect to magnetization 
direction ([100] or [010]) causes a resistance difference 
between these two magnetization orientations. 

While TAMR is a band structure effect, the extreme 
amplification of the effect between 4.2 and 1.7 K can only 
be understood by considering many-particle effects. The 
data of Fig. |2t and b further support the presence of 
electron-electron correlations. 
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FIG. 2: a) Differential conductance-voltage curves at 1.7 K 
including sets of curves belonging to the two magnetization 
states at B=0, and one set at B=300 mT along (p = 3 and 
6°. b)Log-log plot of same data, c) Temperature dependence 
of the zero bias resistance of the high resistance state. 

The differential conductance G(V) = dl/dV in Fig.Eb 
is calculated from current-voltage-characteristics taken 
in the remanent magnetization state after preparation of 
the sample at an angle ip. In Fig. 01, the G(V) curves 
at 1.7 K are segregated into two groups associated with 
each of the two resistance states identified in Fig. [IJi. 
The magnetization of both (Ga,Mn) As layers relaxes into 
the low resistance state corresponding to magnetization 
along [010] for all preparation angles between 45° and 
135° (upper curves in Fig. Et) • It relaxes along the [100]- 
high-rcsistance state for preparation fields within 45° of 
the [100] direction (lower curves). 



The G(V) characteristics have a distinct shape for each 
magnetization direction. High conductance curves are 
square-root-like and show non-vanishing conductance at 
zero voltage. This is typical for the metallic behaviour of 
doped semiconductors (here (Ga,Mn)As) taking part in 
the tunnel process. The square root dependence stems 
from weak localization in metals [T^ . The high resistance 
state, however, shows insulating behaviour. The conduc- 
tance vanishes at zero voltage and the curves follow a 
higher power law as seen in the log-log plot of Fig. \^p. 

The shape of these G(V) curves is reminiscent of those 
observed by Lee et al. [4j in their investigations of tunnel- 
ing from a metal into disordered boron-doped silicon near 
the MIT. Their theory describes the formation of a soft 
ES-gap stemming from an increase in electron-electron 
Coulomb interaction as screening is reduced when the 
MIT is approached from the metallic side. This gap 
around the Fermi energy is observable in the low bias 
tunneling conductance. It manifests itself in a different 
power law behaviour G(V) oc V m , with to = 1/2 for 
metallic (just above the MIT) and to > 2 for insulating 
(below the MIT) material. 

The transition from metallic to insulating behaviour is 
thermally activated as demonstrated by the zero bias dif- 
ferential resistance data taken at an angle corresponding 
to the high resistance state at 1.7 K, shown in Fig. 
The plot shows that the data is consistent with an ex- 
ponential activation following R oc exp (1/T) 1 / 2 expected 
for an ES gap material under the usual assumption of 
single hop tunneling Q, |(| . We note however that be- 
cause the accessible temperature range spans less then 
one decade, and we cannot rule out an R oc exp(l/T) 1 / 4 
dependence expected from a non correlated Mott transi- 
tion. 
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FIG. 3: Resistance of the sample at 1.7 K and 300 mT as a 
function of magnetic field direction under various bias and on 
different cooldowns. d) and e) each show the results of two 
separate measurement, confirming reproducibility. 
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A further clue to the origin of our effect comes from 
the data of Fig. showing the resistance behaviour of 
the sample at 300 mT, a field sufficient to force the mag- 
netization vectors parallel to the field direction for all 
field angles. The figure shows very strong and apparently 
random oscillations reminiscent of quantum interference 
effects, and likely arising from a statistically defined elec- 
tronic state in the sample. Fig.^a-c show that the ampli- 
tude of these oscillations grows significantly with reduced 
bias for angular ranges (gray regions in the figure) asso- 
ciated with ES behaviour. Fig. [3>e show a significant 
effect upon thermal cycling to temperatures of some tens 
of K (above the ES activation temperature). The three 
curves are nominally identical measurements after subse- 
quent thermal cycling. The main behaviour of the device 
remains unchanged, but the details of the fluctuations 
change significantly, presumably corresponding to a new 
impurity configuration. We emphasize that the change 
in the fluctuation pattern is due to thermal cycling. On 
a given cooldown, measurements are reproducible as con- 
firmed by the two, almost indistinguishable, sets of mea- 
surements included in each of Fig.[3|i and e. G(V) curves 
at 300 mT along the <p = 3 and 6° given in Fig. |3; show 
similar behaviour to the other high resistance curves. 

We attribute these fluctuations to quantum interfer- 
ence effects on the variable range hoppin g th at is the 
transport mechanism in the insulating state |l4| . In refer- 
ence jl5| , such quantum interference was confirmed by a 
statistical analysis of the amplitude of the fluctuations as 
the Fermi energy is varied using an external gate bias. A 
similar analysis cannot be reliably applied here as the re- 
lation between the magnetization direction and the Fermi 
energy is complex and non-linear. 

The behaviour of our sample is fully consistent with the 
depleted (Ga,Mn)As injection layer undergoing a MIT 
triggered by a reorientation of the magnetization from 
[010] to [100], and the formation of an ES-gap. In order 
to understand how a change in magnetization direction 
can trigger a MIT, we need to consider the criterion for 
the passage from metallic to insulating properties. 

Metallic transport properties require the charge carri- 
ers to be in long range itinerant states. This condition 
implies significant overlap of the wavefunctions between 
acceptors (Mn atoms). As such, the passage through the 
MIT depends not only on the carrier density, but also on 
the volume occupied by the wavefunctions of the states 
of the individual dopants. Indeed a change in impurity 
conductivity associated with a strain induced change in 
Bohr radius has been reported previously |la. 

For (Ga,Mn)As, it is well established jllLll^ that de- 
tails of the DOS are influenced by the magnetization di- 
rection. To determine if this plays a role in renormalizing 
the wavefunction extent of our localized states, we turn 
to a numerical calculation of bound hole states. 

The valence band in a zinc blende semiconductors is 
described in a 6-band kp framework by the full Luttinger 



Hamiltonian H(k) |rjj. Luttinger and Kohn showed |l8| 
that a charged impurity in a semiconductor host can be 
described by treating the impurity wavefunction as an 
envelope wavefunction B(k) of Bloch waves. The time- 
independent Schrodinger equation is then [Isj : 

(k)B(k) + J d 3 k'U{k - k')B(k') = E b B(k) (1) 

where E\> is the binding energy. U is the Fourier trans- 
form of the Coulomb limit of the Yukawa interaction po- 
tential between the hole and impurity-center, -ff(k) is 
a 6 x 6-matrix acting on the 6 dimensional vector B(k) 
locally in k-space. It describes the bandstructure includ- 
ing, for magnetic (Ga,Mn)As, pd-exchange, strain and 
SO-coupling. Eq.Q is solved numerically for the accep- 
tor ground state in a magnetic zinc blende semiconduc- 
tor. The obtained k-space wavefunctions arc fitted to hy- 
drogen ground state wavefunctions in order to extract an 
effective Bohr radius for each k direction and for various 
magnetization directions. Details of these calculations 
will be published elsewhere. 

From symmetry considerations, one expects that our 
(Ga,Mn)As has four-fold in-plane magnetic easy axes. 
Experimentally however, a small additional [010] uniaxial 
anisotropy is generally observed in samples from various 
sources. While the origin of this uniaxial anisotropy is 
unknown, considerable insight is gained by phenomcno- 
logically introducing a uniaxial strain term e u along [010] 
|8j . An important result of the bulk band structure calcu- 
lation is that the energy ordering of the bands depends 
on the relative strain and magnetization directions. In 
Fig. Et, the T point energy of the top 4 valance bands is 
plotted as a function of the amplitude of e u and the pd- 
exchange. For zero strain, we obtain the expected [lllll2j 
nearly equidistant spacing between the 4 bands, while in- 
creasing strain shifts the relative energy position of the 
light-hole-likc bands with respect to the heavy-holc-likc 
bands, leading to a change in ordering of the bands at 
the line defining the crossing point between the top two 
planes in the figure, and thus to a profound change in the 
character of the wavefunction at a given energy. This is 
directly reflected in the wavefunction of the impurity hole 
state. 

Fig-Eb shows the extent of a bound hole for magnetiza- 
tion along different in-plane easy axes with an amplitude 
of e u beyond the crossing point for the bulk energy at 
the T-point. There is a significant difference in the size 
of the wavefunction between the M||e u and the M _L e u - 
state. The result is that the wave function overlap be- 
tween neighboring dopant sites depends on the direction 
of the magnetization, and for reasonable parameters is 
consistent with the observation of the magnetization re- 
orientation induced MIT. 

It is important to note that the above calculation does 
not predict a MIT for the bulk (Ga,Mn)As which, as typ- 
ical in high quality (Ga,Mn)As, has nearly temperature 
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FIG. 4: a) F point energy of the four top bands as a function 
of strain e u and pd-exchange. b) Extents of a hole bound to 
a Mn-impurity. [001] is the growth direction and e u is along 
[010]. Left: Magnetization (M || [100]). Middle:(M || [010]). 
Right: Perspective view for M || [010]. 

independent resistivity in the temperature range of in- 
terest. Our prediction is limited to the thin (Ga.Mn)As 
layer near the interface with the tunnel barrier. The bar- 
rier consists of LT-GaAs, a material with many mid-gap 
traps and leading to a gradual spatial depletion of the 
(Ga,Mn)As near the barrier on the length scale of the 
Thomas-Fermi screening length of ~ 2A, and thus to a 
much lower effective local carrier density, as illustrated 
schematically in Fig.^. If the thin depleted layer is close 
enough to the MIT, a change in the wavefunction ex- 
tent can trigger the transition. Another point is that the 
mean-field p-d exchange strength is smaller for lower hole 
densities. As seen in Fig.0J a weaker pd-exchange lowers 
the amplitude of e u needed to reach the band crossing to 
much smaller values that are consistent with those used 
in the original description of TAMR HQ. 

A key point here is that while this thin depleted layer 
plays a limited role in the magnetic properties of the 
(Ga,Mn)As, it has a dominant effect on the transport. 
The reason for this lies in the very short mean free path 
of holes in (Ga,Mn)As which is of the order of a few 
A. Combining the fact that, by definition, the transi- 
tion from diffusive to tunneling transport takes place at 
a density where holes can no longer classically diffuse, 
with this extremely short mean free path, must lead to 
a very thin effective injector layer from which the tun- 
neling originates. This injection layer is characterized by 
considerably reduced carrier density relative to the bulk, 
and thus consistent with our model description above. 

The link between the MIT and the development of 
an ES-gap is that the reduction in mobile carriers re- 
duces screening. With unscreened Coulomb-interactions 
between localized states, the DOS at the Fermi-energy 
for carrier concentrations just below the MIT vanishes 
with a dimension-dependent power lawQ, i.e. D(e) oc e, 



where e is the energy measured with respect to the Fermi 
energy, for two-dimensional materials, and D(e) oc e 2 in 
3D. Of course, we cannot rule out that the change in 
screening causes a slight shift in position of the effective 
injector thus modifying the tunneling distance, and ac- 
counting for part of the observed resistance change. 

Larkin and Shklovskii^ derived a power law behaviour 
for conductance versus voltage (G-V) curves for tunnel- 
ing between two three dimensional localized materials 
where the parabolic DOS D(e) oc e 2 on each side mul- 
tiply and lead to G oc V 6 . In our case, given the very 
thin nature of the the injector, a 2-D description of the 
DOS of the tunneling reservoirs is likely more appropri- 
ate. Thus, the DOS is linear in energy, and the expected 
power of the G-V curves is 4. Moreover, it is unclear in 
our experiment whether both sides of the barrier play an 
equivalent role, or whether the effects are dominated by 
one electrode. Depending on the relative role of the two 
barriers, one would thus expect a power law somewhere 
between oc V 2 and oc V 6 . Neglecting very low voltages 
where thermal smearing is important, this prediction is 
in good agreement with the experimental observation of 
Fig. El 

In conclusion, we have observed the magnetization re- 
orientation induced MIT in a (Ga,Mn)As based transport 
sample. The transition can be understood as stemming 
from a modification of the wave functions of individual 
dopants due to the coupling of the Mn dopants to the 
magnetization direction in the bulk. Furthermore, the 
transition is accompanied by the opening of an ES-gap 
at the Fermi energy which manifests itself as a change in 
the power law behaviour in conductance-voltage charac- 
teristics in tunneling experiments. This is the first obser- 
vation of a MIT induced by a change in magnetization 
direction in any material. In addition to the fundamental 
interest intrinsic to these observations, the results may 
also have technological relevance in opening up new pos- 
sibilities of controlling the transport properties of devices 
by magnetization reversal. 
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